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1 Introduction 

Construction of the massive gravity is plagued with theoretical problems and issues 
that have not been resolved from the pioneered work of Fierz and Pauli |lj q One 
such a well known problem of massive gravity is that there is no smooth massless 
limit in the perturbation theory of massive gravity in the sense that the massless 
limit in massive gravity exists but does not agree with the results derived from Gen- 
eral Relativity (GR) that describes massless gravitons. This pathological behavior 
of massive gravity is known as van Dam-Veltman-Zakharov discontinuity [6l [7]. 

Even if the construction of massive gravity is interesting theoretical challenge 
there is another more stronger reason for their formulation. In fact, generalizations 
of GR with a small, non-zero, graviton mass leads to large scale (or infrared) mod- 
ifications of General Relativity. Since GR has been directly tested from scales of a 
fraction of millimeters up to Solar System scales it is possible that the structure of 
the theory changes on the large distances and hence infrared deviations from GR 
cannot be excluded. Then the massive gravities could modify GR at large cosmolog- 
ical scales and also can explain recent accelerated expansion of the universe without 
assuming the existence of mysterious dark matter and dark energy. 

At present it is not clear how to formulate consistently a theory of massive 
gravity. Since the Einstein's theory of gravity is well tested theory it seems to be 
natural to add to the action of GR a term which will, in the linearized approximation, 
give a mass to gravitons without modifying the kinetic terms coming from GR. 

Very interesting formulation of the massive theory of gravity is based on Brout- 
Englert-Higgs mechanism for gravity. In more details, we consider a diffeomorphism 
invariant action with usual Einstein-Hilbert term together with the function of the 
metric coupled to D scalar fields for D— dimensional gravity [8]. Then gravitons 
acquire a mass due to a mechanism which may be thought of as the Brout-Englert- 
Higgs mechanism for gravity when the vacuum expectation value of each scalar field 
breaks one coordinate reparameterization invariance q 

Our goal is to develop the Hamiltonian description of the systems studied in 
m m [ini [m Ha [la [H]. with the help of (D - l) -l- l spUt formalism for the 
space-time metric we find corresponding Hamiltonian and show that it is linear 
combination of the Hamiltonian and diffeomorphism constraints. The consistency 
of the theory demands that these constraints are preserved during the time evolution 
of the system. In order to check whether they are preserved or not we calculate the 
Poisson brackets of these constraints. However it turns out that it is non-trivial task 
to calculate the Poisson bracket of the Hamiltonian constraint for the scalar field 
whose dynamics is governed by general action. Despite of this fact we show that 
the Poisson brackets of the scalar fields Hamiltonian constraints are proportional to 
the diffeomorphism constraints. In other words we show that the Poisson bracket of 
the scalar field Hamiltonian constraints takes exactly the same form as the Poisson 



^For extensive reviews of various aspects of IR modification of gravities, see [5] |31 IH [S] 
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bracket of General Relativity Hamiltonian constraints. Collecting all these results 
we obtain the Poisson algebra of constraints is closed and hence these constraints 
are consistent with the time-evolution of the system. This result is crucial for the 
possibility of the fixing the gauge in the framework of the Hamiltonian formalism. 
We fix the gauge by introducing D gauge fixed functions that correspond to the 
gauge fixing conditions introduced in [U [U [TOl [HI [121 [131 [II]- Then the system of 
these gauge fixing functions and original constraints form the set of the second class 
constraints that can be explicitly solved. Further, since the original Hamiltonian 
is given as the linear combination of constraints we find that in the process of the 
gauge fixing it strongly vanishes. On the other hand the gauge fixing implies that 
the reduced phase space is spanned by the spatial components of the metric hij and 
their conjugate momenta p^K Then we argue that the the gauge fixing condition 
t = (f)^ naturally introduces the Hamiltonian on the reduced phase-space equal to 
— Po where po is the momentum conjugate to (fP . Note that po is the function of 
reduced phase space variables h^j , p*-' as a result of the solving of the Hamiltonian 
constraint. Say differently, we claim that the definition of the massive gravity is 
given by the reduced phase space variables hij,p^^ together with the gauge fixed 
Hamiltonian Hfix = —po{hij,p^^), and where all non-dynamical modes are absent. 

Experiences from many areas of theoretical physics teach us that in some cases 
the Lagrangian formulation of given theory is much more efficient then the Hamil- 
tonian ones. For that reason we feel that it is useful to determine the Lagrangian 
formulation of the gauge fixed theory as well. It turns out that in order to find 
this Lagrangian it is convenient to introduce new non-dynamical modes so that the 
Legendre transformation from the Hamiltonian to Lagrangian formulation can be 
easily performed. Interestingly this Lagrangian can be written in such a form that 
resembles the standard Einstein-Hilbert action together with the potential term that 
explicitly breaks the diffeomorphism invariance of given theory. 

The structure of this paper is as follows. In the next section (jSj) we perform the 
Hamiltonian analysis of the system introduced in [H [9l [TOl [HI [121 IISI HI] ■ Then 
in section ([3|) we calculate of the algebra of constraints. Section (jlj) is devoted to 
the study of the gauge fixed theory. In section ([5]) we give two examples of scalar 
potentials that allow us to find explicit form of the gauge fixed Hamiltonian. Finally 
in conclusion dHj) we outline our results and suggest possible extension of this work. 

2 Hamiltonian Analysis of General Relativity with 
Massless Scalars 

In this section we develop the Hamiltonian formalism for the following action 

^ = 16^ / ^""^v^t^ - ^(^^'')] ' (1) 

where G is D— dimensional Newton's constant and the induced internal metric is 
defined as 

H^B ^ g,u^^^A^^^B ^ (2) 



where (j)^ are real D scalar fields with A = 0, . . . ,D — 1 and where the indices 
A,B,... are raised and lowered using the metric tjab = diag(— 1, 1, . . . , 1). Finally, 
L is apriori a generic function of H^^. The detailed analysis of properties of function 
L was performed in [101 [HI [131 [II]- It is important to stress that the form of the 
potential L is not arbitrary. In fact, it was argued in [TUl JUl [131 [E] that this 
potential has to lead to the equations of motion that possess following vacuum 
solution 

0^ = x>'S^ , g^, = T]^, . (3) 

This requirement leads to the constraint on the potential L 

^HaT = 2VabL{H,) , (4) 

where H^^ = rj^^ . 

Our goal is to develop the Hamiltonian formalism for system defined by the action 
([T|). As usual in the study of the Hamiltonian formalism for gravity we introduce 
(D — 1) + 1 formalism. Explicitly, let us consider D dimensional manifold M. with 
the coordinates x^ , /i = 0, . . . , D — 1 and where x^ = (t, x) , x = (x^, . . . , x^~^). 
We presume that this space-time is endowed with the metric g^y{x'') with signature 
(— , +,...,+). Suppose that M. can be foliated by a family of space-like surfaces 
Sj defined by t = x°. Let /ijj, i,j = l,...,D denotes the metric on S^ with inverse 
h^^ so that hijh^^ = 6f. We further introduce the operator Vj that is covariant 
derivative defined with the metric hij. We define the lapse function N = 1/ ^J—g^^ 
and the shift function A^* = —g^^/g^^- In terms of these variables we write the 
components of the metric g^y as 



9m = 
„oo _ 


1 


+ Nih'^Nj 




9oi 

1 


9'' = 




9ij 




g 


m ' 


to see 


that 

















(5) 



^ = iVVdet h . (6) 

In the {D — 1) + 1 formalism H"^^ takes the form 

H^"" = -V„0^V„0^ + h^W,(l)%<P^ , (7) 

where 

V„0^ = ^W^-Ar^9,0^). (8) 



Then from ([T|) we easily find the momenta conjugate to 



kA 



Vdeth 5L ^ ,r. 



Using this result we find following matrix equation 

Kab = Dac{H^''W'' - H^'^)DnBiH^'') , 

(10) 

where we introduced following matrices 



5Hab 

(11) 

For further purposes we also introduce the matrix D^^ inverse to Dab so that 
DabD^"^ = S^. 

Now we presume that the equation fITU]) can be solved for Hab- Let us denote 
this solution as J"^^ = J'^^i^K, V). Then (TTOj) implies following relation 

^yAB _ jAB^ ^ D^''{J)KcdD''''{J) (12) 

that will be useful below. Collecting all these results we find the Hamiltonian of the 
scalar fields in the form 

H'^= I d''i^{dt<p^pA -L)= I d''^{N{^)nU^) + N\^)nf{^)) , (13) 



where 



n^ = -^(2Vd^KABD''^ + V{J) 



16nG 
nf = PAdi< 

(14) 



^«r ) 



and where d = D — 1. 

In the same way we should proceed with the Hamiltonian analysis of the General 
Relativity action. Since the procedure is well known we immediately write the final 
result 



H^"^ = I t/'^x(Ar(x)7{^^(x) + ArXx)-Hf ^(x)) , (15) 



where 



^r = ^^^'^Q^,u^'' - ^|^i?(^-i) , n, = -2K,V,7:^' . (16) 

Vdet h IbvrG 

Let us explain notation used here, tt*-' are momenta conjugate to hij. The generalized 
metric Q"^^^^ is defined as 

Qijki = T^ihikhji + hiihjk) — hijhki . (17) 



The inverse metric Q^^^'' is equal to 

gijkl ^ \^^^kyl ^ ^lyk-^ _ ^ij^kl _ ^^g^ 

Finally, w-^~^\ Vj are Ricci curvature and covariant derivative calculated using the 
metric hij. Finally due to the fact that the action ([1]) does not contain time deriva- 
tive of A^, A^* we find that corresponding conjugate momenta tttv, vtj are primary 
constraints of the theory 

-Ki^Q , vr^ ^ . (19) 

The condition of preservation of these constraints implies the existence of the sec- 
ondary ones 

Ht = n^^ + H^ ^ , n, = nf^ + nf^o (20) 

or their smeared form 

Ht{M) = I d'^^M{^)'HT{^) , Hs{M') = f d'^xM'{x)n^{^) , (21) 

where M(x), M*(x) are arbitrary functions. 

3 Algebra of Hamiltonian Constraints 

As the next step we have to demonstrate the stability of the constraints fl20|) . or 
equivalently, we have to show that these constraints are preserved during the time 
evolution of the system. In fact, the careful analysis of these constraints was per- 
formed in 70's in the geometrodynamics program fiEi [T6l [T71 [18] . It was argued 
there that for the consistency of theory the Poisson brackets of constraints should 
have the form 

{Hs{N'),Hs{M')} = Hs{N'd,M^ - M'd.N^) , 
{HsiN'),HTiM)} = HriN^d^M), 
{Ht{N),Ht{M)} = HsiiMdjN - NdjM)h^') . 

(22) 

The aim of this section is to show that the Poisson brackets of constraints fl20l) obey 
the algebra fl2^ . In fact, it is well known that the General Relativity constraints 

H^^{N) = f d^xN{x.)n^^ix) , H^^{N')= f d^icN'{^)nf^ix) (23) 

obey the Poisson brackets relations fl22|) . Further in case of scalar field diffeomor- 
phism constraint we easily find 



H^{N'),Hl{M') \ = H^{N'diM^ - M'diN^) . 

(24) 



Then using the fact that the mixed Poisson brackets < Hg{N^), H^^{M^) > are triv- 
ially zero we find that the Poisson brackets of diffeomorphism constraints obey the 
first equation in ( l22l) . 

As the next step we calculate the Poisson bracket between diffeomorphism con- 
straint Hs{N^) and Hg(M). Using following Poisson brackets 

{Hs{N^),Pa} = -d^iN'PA), 

{Hs{N'),h,^} = -N''dkK,-d,N''hk,-Kkd,N\ 

lHs{N'),^/Aeih\ = -N'^dkVdet h - dkNWdet h 

(25) 

we easily find 

{HsiN'), Kab] = -N'^dkKAB , {HsiN'), V^""} = -N'^dkV^'' . (26) 

Collecting all these results we obtain 

\^HsiN'),n^^ = -N'^dk'Hi^ - dkN'^n^ . (27) 

Finally using the fact that {H^^{N'), H^^{M)} = H^^{N'diM) and the equation 
fl27j) we find that 

{HsiN'),HTiM)} = HriN^d^M) . (28) 

As the last step we calculate the Poisson brackets of the Hamiltonian constraints. 
While it is well known that the Poisson bracket of the GR Hamiltonian constraints 
takes precisely the form given in (l22l) it is far from obvious that for the system defined 
by the action ([1]) it holds as well. In case of simple scalar action ~ g^^d^^'dycf) the 
result is well known, see for example [Ul EHl [IZl [IE] ■ Further, it was shown very 
elegantly in [TJ] that the scalar field action in the form L{—g'^'^d^(j)di,(j)) leads to 
the Hamiltonian formulation where the constraints obey the Poisson brackets (!22|) . 
On the other hand the scalar field Hamiltonian constraint ([H]) is more general then 
the form of the Hamiltonian constraints studied in [15] so that we perform the 
explicit calculation of the Poisson brackets of the smeared form of the Hamiltonian 
constraints (TT^ below. 

Before we proceed to this calculation we discuss the calculation of the mixed 
Poisson brackets between the GR Hamiltonian constraint and scalar fields Hamil- 
tonian constraint. The crucial point is that the constraint Hj, depends on gij only. 
Then it is easy to see that 



H^^{N),H^{M)\ + Ih^{N),H^^{M)\ = 
I d^^d^yN{^)M{y) ({?/^^(x),?/^(y)} + |?/^(x),?/?^(y)} 



using the fact that |?{^(x),'H^^(y)| ~ (5(x - y). 



(29) 



3.1 Poisson Brackets of Hamiltonian Constraints for Scalar 
Field 

In this section we determine the Poisson bracket between Hj,{N),Hj,{M). As a 
warm example we begin with the single scalar field action 



S = d^Xy/^L{-I/2) , (30) 

where / = g^'^d^(j)dp(j). In (Z^ — 1) + 1 formalism this action takes the form 

S= j d^xNVd^L{-{Vl(p - h'^dicpdj^)) . (31) 

Then it is easy to see that the momentum conjugate to (f) is equal to 



p = Vdet hL'{I/2)Vn<P (32) 

and consequently 



K 



P 



Vdet h 



L'\I/2)iI + V), (33) 



dLix) 



where V = h'^Wi(f)dj(j) and where L'{x) = ^^ . Let us presume that the equation 
( 133|) can be solved for I and we denote this solution as / = J{K,V). Then (l32l) 
imphes 

Vn0 = -^ TT (34) 

and hence the Hamiltonian takes the form 

H= I rf'^x(Ar(x)HT(x) + N\^)'Hi{^)) , (35) 

where 

p2 



Ut = Vdet hL(J/2) . (36) 

Vd^L'{J/2) 

Our goal is to calculate the Poisson bracket < Hj,(N), Hj^{M) >. To do this we start 
with the calculation of the following Poisson bracket 

2p 1 



{4>.Ht] = Jd''^{(j),N{^)HT{^)} = N 



^A^L'\J/2) 



where we used 



- N P' L"iJ/2) 5J p 5J 



^^(-)'^^^)>-2dl7ll^^-)^^--^) 



(37) 



(38) 



To proceed further we note that J obeys the equation 

K = L'\JI2){J + V) . (39) 

When we difference this equation with respect to K and use the fact that J = 
J{K, V) we find 



KL"iJ/2)^ = L'iJ/2) (l - L"{J/2)^ \ 



6J\ 

(40) 

Inserting this result into fl37j) and performing the appropriate manipulation we find 
that the Poisson bracket fl371) takes the form 



(41) 



{0, HriN)} = N ^ ■"" = ATV, 



that agrees with the result derived in [T5] . 

As the next step we determine the Poisson bracket between p and Ht{N) 



{p,Ht{N)} = I d'^^N{x) X 



^' -L"(J/2)^(x) {p, \/(x)} - lVd^L'iJ/2)^M {P, ^(x)} 



2Vd^L'^{J/2) dV '^-^ ^ '^ 2 'dV 



With the help of the relation f p9|) we find 

dJ T'3fj/n\ A , dJ\ dJ T,dJ 



(42) 



KL'V/2)^--L'V/2){l^^). ^--L'jj^- (43) 

Inserting these expressions into (H2|) we find 

{p,i7T(A^)} = - [ d'^xNVdethL'ix.) {p,V{x.)} = d^[N ^ det hh'^ L'd. 



JY\ ) 



(44) 

where we also used 

{j9(x), V{y)} = -20{y)dy,<P{y)dy.6{^ - y) . (45) 

Now we are ready to determine the Poisson bracket {Ht{N),Ht{M)}. To do this 
we consider following expression 

{{0, Ht{N)} , Ht{M)} - {{0, Ht{M)} , Ht{M)} . (46) 



In the first step we use (HTj) and we find 

N , ■"" ,Ht{M)\- \m ^ ■"" ,Ht{N)\ . (47) 

Clearly all terms in f H7|) that are not proportional to the partial derivatives of A^, M 
cancel. To proceed further in the calculation of (H7|) we have to calculate following 
Poisson bracket 

X (-L'(J/2)h'^d,4>dA , ^ ] + -^di[Vd^h'^di(t)]] . 

V '^Vd^L'{j/2y deth '^ '^^ J 

(48) 

Then it is easy to see that ^7^ {M{^,Ht{N)} -N{^,Ht{M)}) = 0. Finally 
with the help of (l44l) we find that ( 146!) is equal to 

{{0, Ht{N)} , Ht{M)} - {{0, Ht{M)} , Ht{N)} = {Nd,M - MdiN)h'^d,<j) . 

(49) 

With the help of the Jacobi identity we can rewrite this expression into the form 

{(J),{Ht{N),Ht{M)}} = {NdiM -MdiN)0dj(f) = 

= {(t>,Hs{Nd^M-Md^h'') + f{<P)] 

(50) 

for arbitrary function /(0). On the other hand performing the same step with vr 
we find that / does not depend on cf) as well and hence could be taken to vanish 
at least at classical level. In summary, we proved that the Poisson bracket of the 
Hamiltonian constraints of the scalar field with general action fl30|) takes the desired 
form 

{Ht{N), Ht{M)} = HsiiNdiM - MdiN)h'^) . (51) 

Now we are ready to proceed to the calculation of the Poisson bracket of the Hamil- 
tonian constraint when the action for the scalar fields is given in ([T]). We again start 
with the calculation of the Poisson bracket between (j)^ and Hj, 

oKmn oKmn J 

(52) 

To proceed we note that the equation (iT2ll implies following relation 

Kab = Dac{V^'' - J'''')Ddb , (53) 



where Dab depends on J. Then we calculate the Poisson bracket between (j)^ and 
the relation given above, multiply the result with with D^^ and finally take the 
traces over capital indices. As a result we find 



D^^ {0^, Kab} = ^D^^j^ {0'', Kmn} D^^'K^b + i^DnB {0^, Kmn] • 

oKmn oKmn 

(54) 
Using this result in (152|) we find 



(55) 



In the similar way we proceed with the calculation of the Poisson bracket between 
Px and Hj, 



(56) 
where we used the Poisson bracket between px and ( l53l) that implies 



(57) 
As in case of single scalar field we consider following expression 

{{<f>^H^{N)} hUm)} - {{<f>\H^{M)} H^N)} , (58) 

where the first double Poisson bracket is equal to 

[[<P^,H^{N)}hUm)} = -^D^^d,[MVd^h^WBcdj<P''] - 

(59) 

To proceed further we have to determine the Poisson bracket between Dab and 
Hj,{N). Note that when we compare the variation of ( l53l) with respect to Kmn 
with the variation of (153!) with respect to Vmn we find following relation 

§jMN 5J^'^^ 

Then using this relation in fl59|) and after some manipulation we find 



||0^,/7^(Ar)} /7^(M)} - ||0^,i7^(M)} i/^(Ar)} = (AT^iM - MdiN)0dj(j)^ . 

(61) 
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Performing the same analysis with conjugate momenta pa we find the desired result 

{h^{N), H^{M)\ = H^{{NdiM - MdiN)0) . (62) 

The upshot of this long analysis is the proof that the Poisson brackets of Hamiltonian 
constraints Ht{N) = Hj!^{N) + H^{N) take exactly the same form as in fl22|) . 
Using this fact we immediately obtain that the Hamiltonian and diffeomorphism 
constraints are preserved during the time evolution of the system. This is very 
important result since only after determining the complete constraint structure of 
given theory it is possible to perform the Hamiltonian gauge fixing. 

4 Fixing Gauge 

In this section we fix the gauge freedom that in the Hamiltonian treatment are 
expressed by an existence of D first class constraints 'Ht{'x.) ~ 0,'Hj(x) ^ 0. This 
procedure is an analogue of the Higgs mechanism used in the construction of the 
massive gravity O [101 [HI Ull [13] . In these models the vacuum expectation values of 
the scalar fields coincide with D space-time coordinates. The result of this fixing is 
the complete breaking of the space-time diffeomorphism and emergence of the mass 
term for the graviton in the action when we study the small fluctuations of gravity 
above the flat space-time. 

The standard way how to flx the gauge freedom in Hamiltonian framework is to 
introduce D gauge flxed functions that can be interpreted as additional constraints 
imposed on the system and that have non-zero Poisson brackets with the original flrst 
class constraints. As a result the extended system of original constraints together 
with gauge flxed functions form the collection of the second class constraints with 
no gauge freedom left (For review of this formalism, see [191 1201 [21].). 

With analogy with the flxing the gauge given in [SI [IDl [HI [121 [13] we introduce 
following D gauge flxing functions 



6^^(x) = </>^(x)-a;^. (63) 



Clearly 
and 



Due to the fact that the Poisson brackets (IMI) and (1651) are non-zero on the constraint 
surface we see that the collection of constraints {Q^,T-LT,T-ii) is the system of the 
second class constraints. Alternatively, the requirement of the time preservation of 
the constraints ^"^(x) ^ during the time evolution of the system implies following 



11 



consistency equation 



dt - ^^^W'^^ 



^ ^ ^ ^ 2deth 6KiN(x.y ^ ' 



"^^'^""^ = d,g\^) + {g\^),H} 



dt 



zdetn. dKQN\y) 



(66) 



where we used ( IMl) and ( l65l) . We see that these equations determine A^ and A^* as 
functions of the canonical variables hij,p^^ with no gauge freedom left. 

The fact that Q"^, T-Lt, "Hi are the second class constrains implies that they vanish 
strongly and can be explicit solved for pa- As a result the reduced phase space is 
spanned by hij and p'^K Further, since the Hamiltonian of the original system was 
given as a linear combination of the constraints Ht, 'Hi we now see that it vanishes 
strongly. 

On the other hand let us write the original action ([T]) in the form 



S= dtd''^{dth,'K^' + dt(t>''pA-H) = dt{ d^^idtK.-K^'+p^ip'^g.i))) , (67) 

where we used the fact that H = and imposed the gauge fixing functions (!63|) . 
We see from (!67|) that it is natural to interpret —po as the Hamiltonian density of 
the reduced theory 

^/.. = -po(/^..(x),H?^(x),Hf«(x)) ^ -^^^Po , (68) 

where we also used the fact that from the constraints Hj = we can express pi as 

p.(x) = -^^(x) . (69) 

Finally note that po can be derived from the Hamiltonian constraint Tir = 'Hj.^ + 
'H|.(x) = at least in principle. 

In summary, we found the Hamiltonian formulation of massive gravity where the 
physical degrees of freedom are h^j , p*-' and where the Hamiltonian is given in f lBS]) . 
Note also that the explicit form of this Hamiltonian depends on the form of the 
function L[H^^). We give simple examples of two solvable potentials in the next 
section. Generally however it is very difficult to find explicit form of the gauge fixed 
Hamiltonian due to the complicated structure of the function L{H^^). 

Despite of this fact we now show that it is possible to find the Lagrangian density 
for given gauge fixed theory. To to this we introduce four modes A, B, Ci, D^ and 
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corresponding conjugate momenta {pA_,PB,p\Pi) with non-zero Poisson brackets 

{A(x),p^(y)} = ^^(x-y), {5(x),pB(y)} = 5(x-y), 
{a(x),p'(y)} = 6i6{^-y), {D\^),p,{y)} = 6]6{^- y) . 

(70) 

With the help of these additional modes we rewrite the Hamiltonian for gauge fixed 
theory as 

^/« = -poiK^,Aa)+B{n^''-A)+D\nf''-a)+VAPA+VBPB+v%+v,f , (71) 

where the Lagrange multipliers va,vb,v\ Vi ensure that pa, PB,Pi and j9* are primary 
constraints of the theory 

Pa^O , Pb ~ , Pi^O , p' ^0 . (72) 

Then the fact that these constraints have to be preserved during the time evolution 
of the system implies the secondary constraints 

dtPA = {pA,Hf,,} = ^-B = ^A-0, 

dtPB = {pA,Hfi,} = -A + n^'^ = <^B^O, 



dtp' = {p\Hf,,} = ^-D' = <l>'^0. 



(73) 



It can be shown that the collections of the constraints {pA,PB,Pi,p\ ^a, ^b, ^\ ^i) 
are the second class constraints. Solution of these constraints reduces (17T|) into the 
original form of the gauge fixed Hamiltonian (168!) . 

The main advantage of the extended form of the Hamiltonian density 07ip is 
that it allows us to find corresponding Lagrangian in relatively straightforward way. 
In fact, from (I7T|) we easily obtain the time derivatives of the canonical variables 
hij, A, Jd , U , Cj 

or) /^ 

V det h 
dtA = {A, Hfi^} = va , dtB = {B, Hfi^} = vb , 



(74) 



It turns out that it is useful to introduce following object 



Ki, = ^{dth,, - V,D, - V,A) (75) 
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that due to the first equation in (^^ is related to vr*-' as 

IGtiG „ ^i 



Kij = — - GijkiT^ ■ (76) 

V det h 



Then it is easy to find corresponding Lagrangian 



^f^. = I d'^l T^^^^^^'^^'"'^'^' + ^^"""'^^ + ^^^Poih^,' ^' C,) + BA + D^C, 




' 167rG SttG 



where g is D— dimensional metric with components 

^00 = -5' + D,WW^ , ^0^ = ^^ , 91, = hij . (78) 

and where Ris D + 1 Ricci scalar built from this metric. We see that the last form of 
the Lagrangian (177j) can be interpreted as the sum of the General Relativity action 
with additional potential terms that breaks the full diffeomorphism invariance of the 
theory. It is important to stress that this potential term depends on the auxiliary 
fields A and Cj. In principle these terms could be integrated out however we expect 
that the resulting Lagrangian would be very complicated. 

5 Examples of Potentials L{H^^) 

In this section we give two solvable examples of the scalar function L{H^^) that 
allow to find the explicit form of Lagrangian for massive gravity. 

In the first case we follow [8] and consider function L{H^^) in the form 

L = A + Habv''^ ■ (79) 

Then it is easy to see that 

SL , , 

TTT- = "^BA (80) 

and 

J^^'VBA = V^^VBA - K^^'VBA ■ (81) 

Using these results it is straightforward exercise to find scalar field Hamiltonian 
density 

^T = ^^ (Kabv''^ + V^^VBA + A) . (82) 

iovrG 

Then following the general procedure outlined in previous section we find 



2 
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(77) 



Finally we find the Lagrangian density of the gauge fixed theory in the form 



IQttG 

(84) 

We would like to stress that we can integrate out auxiliary fields A,B,C^,Di from 
the Lagrangian fl8^ and then to derive the Lagrangian density for dynamical modes 
hij only. However the resulting Lagrangian would be very complicated and hence 
we prefer to work with the extended Lagrangian f l8^ . 

As the second example of exactly solvable theory we consider the Lagrangian 
function L{H^^) in the form 

L{Hab) = a + Vfi + H^'^VBA , (85) 

where A and Q are constants. For (l85l) we easily find 

Dab = 1^=^§== , D^^ = 2r^^^V^ + H^B^sA (86) 

2 y/fi + H^^r]BA 



-^ VBA = 1 , ,r^ nA • i^') 



and hence 

tAb^ _ V^' ^Vba - ^^Ka bv''^ 

1 + AKabv' 

After some calculation we obtain the Hamiltonian density for the scalar field in the 
form 

«* = ^ Vfi + I'-'teWl + 4/aB.,- - #pA . (88) 

IottG IottG 

Following the analysis presented in previous section @ we find the gauge fixed 

Hamiltonian density in the form 






(89) 
and corresponding Lagrangian 



I^/ix = j d''^(-^-^^-deig{R + A) + BA + D'Cr 



Vd^l /167rG\' ,..^ / 167rG \ 



n + h^^5ji 



(90) 



In this section we gave two explicit examples of Hamiltonians and Lagrangians 
for massive gravity. Clearly it would be desirable to understand properties of these 
models further. 
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6 Conclusion 

This paper was devoted to the study of the Higgs mechanism for gravity from the 
point of view of the Hamiltonian formahsm. We performed the fixing of the space- 
time diffeomorphism and we argued that the resulting Hamiltonian corresponds to 
the Hamiltonian of the massive gravity. The my advantage of our approach is that 
this theory is defined on the reduced phase space spanned by the physical degrees 
of freedom hij,p^^ only. On the other hand the price we pay for this property is that 
it is difficult to find the form of the gauge fixed Hamiltonian for general potential 
L{H^^). In fact, we are not able to find explicit form of the gauge fixed Hamiltonian 
for the specific form of the scalar actions introduced in [9l [TOl [HI [121 [131 [H] • On 
the other hand introducing additional auxiliary fields we can determine Lagrangian 
for massive gravity that has the form of the ordinary General Relativity action with 
specific potential terms that break diffeomorphism invariance. Then we can ask the 
question how this Lagrangian is related to the original Lagrangian where we fix the 
gauge as in [9l HOl [Til [I2l HSl Hi] . One can hope that these actions could be related 
by some fields redefinitions. However finding this redefinition seems to be very 
complicated due to the presence of the auxiliary fields A, Ci in the Lagrangian (1771) 
whose explicit integration out would lead to very obscure form of the Lagrangian. 

The next important step in our investigation would be to analyze the spectrum of 
fiuctuations around the fiat space-time background. It would be also very interesting 
to study the classical solutions corresponding to these form of massive gravities. We 
hope to return to the analysis of these problems in future. 
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